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ABSTRACT 

We present initial results and source lists of variable sources in the Large Magellanic Cloud 
(LMC) for which we detect thermal infrared variability from the SAGE (Surveying the Agents of 
a Galaxy's Evolution) survey, which had 2 epochs of photometry separated by three months. The 
SAGE survey mapped a 7° x 7° region of the LMC using the IRAC and the MIPS instruments 
on board Spitzer. Variable sources are identified using a combination of the IRAC 3.6, 4.5, 5.8, 
8.0 /im bands and the MIPS 24 fim bands. An error-weighted flux difference between the two 
epochs is used to assess the variability. Of the ~ 3 million sources detected at both epochs 
we find ~ 2, 000 variable sources for which we provide electronic catalogs. Most of the variable 
sources can be classified as asymptotic giant branch (AGB) stars. A large fraction (> 66%) of the 
extreme AGB stars are variable and only smaller fractions of carbon- rich (6.1%) and oxygen-rich 
(2.0%) stars are detected as variable. We also detect a population of variable young stellar object 
candidates. 



Subject headings: stars: variables: other — stars: AGB and post-AGB 
formation — infrared: stars — galaxies: individual(LMC) 



stars: mass loss - stars: 



1. Introduction 

The Spitzer survey. Surveying the Agents of a 
Galaxy's Evolution (SAGE) of the Large Magel- 
lanic Cloud (LMC) provides an unprecedented op- 
portunity to detect thermal infrared variability of 
the infrared stellar population of the LMC. Opti- 
cal variability studi es of the LMC by the MACHO 
( Alcock et al]|l996l) and OGLE (jPaczvnski et al 



variable s tars in the LMC, e.g. the long period 



19941 ) monitoring projects have revealed the pe- 
riod and luminosity relations for a wide variety of 
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variables (jPraser et al.l 120051 ) . The three month 



time span between the two epochs of observations 
for the SAGE survey can detect and constrain the 
variability of long period variables in the evolved 
star population and young stellar object candi- 
dates. Thermal infrared variability of such objects 
has been studied in the Galaxy with IRAS and 
ISO, and Spitzer's improved sensitivity enables 
such studies in nearby galaxies. The 10 month 
IRAS mission surveyed most portions of the sky at 
least twice a year, separated by months and thus 
also constrained variability of the brighter infrared 
sources in a similar way that the SAGE data sam- 
ples the LMC's stellar populations. 

Most of the sources detected as variables in the 
thermal infrared by IRAS or I SO are evolved stars 
i n our Galaxy. Por example. iHarmon fc Gilmord 
(|l988l) used the IRAS VAPQ index to determine 



^The IRAS VAR index is a probability that the source is 
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the nature of the IRAS stellar population that de- 
lineates the Galactic Bulge. In a statistical way, 
they simulated the periods that the IRAS data 
were sensitive to and found the data were most 
sensitive to long period variables, P> 400 days 
up to 1,400 days. Essentially, large amphtude 
and long period variations are easier to detect 
and are also associated with larger mass-loss rates. 
They used this period to estimate the age and 
mass range for this evolv ed stellar population in 
the bulge. Later work by van der Veen fc Habind 
(1990) also used the IRAS VAR to select and study 
OH/IR stars for ground based measurements of in- 
frared light curves and find that these IRAS vari- 
ables fit well into the period-luminosity relations 
found for Mira variables in Baade windows, con- 
firming that OH/IR stars are more evolved forms 
of Miras. ISOGal-DENIS studies of the Galac- 
tic bulge included one field with mul tiple ISOCam 



imag es separated by 17-23 months (jOmont et al 
19991 ). This study revealed a population of red gi- 



ants with weak mass loss for which the brightest, 
and most dusty red giants have variation at 7 and 
15 /xm. These sources are "intermediate aymptotic 
giant branch (AGB) stars" i.e. between early AGB 
and thermally pulsing AGBs, and have low mass 
loss rates, 10'^° to 10"" Mq yr"^ 

In the LMC, the evolved stars are too faint 
for IRAS and ISO measurements to reliably de- 
tect variability. However, IRAS- and ISO- de- 
tected evolved stars have been followed up by 
ground based nc ar-IR var i ability monitoring pro- 



grams. Whitcl ock et al.l (j2003f ) monitored the 



IRAS-selected sample of AGB stars studied by 
ISO and find that larger amplitude variations in K 
are associated with redder stars. The IRAS and 
ISO photometry are sensitive measures of dusty 
mass-loss of the AGB stars. Luminosity is corre- 
lated with period and longer periods have higher 
amplitude variations. Thus, the thermal infrared 
variability of the evolved stars can cause a system- 
atic scatter in the mass-loss rate determinations of 
these evolved stars in the LMC and in mass-loss 



rate vs. luminosity relations (jGroenewegen et al 
20071 ). 

In contrast to the evolved star populations. 



variable in the IRAS 12 and 25 micron photometry, e.g. 
VAR > 50 means larger than 50% chance that the source 
is variable. The variations have to be correlated, i.e. in the 
same direction. 



there have been few thermal infrared variability 
measurements of young stellar objects (YSOs). 
IRAS measurements of variability of YSOs were 
at best difficult due to the low spatial resolu- 
tion of IRAS and the confusion with diffuse ISM 
emission. Nevertheless, for isolated YSOs there 
was a detection of variability at 12 and 25 /im 
for t he Herbig Ae/Be stars AB Aur and WW 
Vul (jPrusti fc MitskevichI \l994 . Intermediate 
mass pre-main sequence stars like UX Ori stars 
show l arge photometric and polarimetric vari- 



ations (jNatta et al.l 119971 ). A concurrent ISO 



and ground based optical monitoring program of 
UX Ori type stars, S V Cep, was perform ed by 



Juhasz et al.l (|2007[ ). lJuhasz et al.l (|20071) found 



that the mid-infrared flux variations were anti- 
correlated with the optical variations, but that 
the far-infrared flux variations were positively cor- 
related with the optical variability. They used 
the infrared variability to help discriminate be- 
tween disk or disk-|-envelope models and find 
that disk models are in b e tter a greement with 



the data. [Abraham et all (|2004l ) used ISO to 



monitor the long term IR evolution of seven FU 
Orionis stars. They detected variation in four 
sources, which provides ten tative support for the 
Hartmann fc KenvonI ( 19851 ) models for accretion 



outbursts in FU Ori objects as opposed to pre- 



cessing jets ( Herbig et al. 20031 ). Interestingly, 



more luminous stars which appear to have similar 
characteristics as the FU Orionis stars have no 
detected outbursts. Chromospheric hot-spots or 
accretion disk variability can cause optical vari- 
ability in pre-main sequence stars like T Tauri 
stars, but these types of objects have very low 
infrared excesses and are unlikely to be in our list. 
Our list would include earlier stages of YSOs that 
have more substantial infrared excesses. Thus, 
while difficult to detect, thermal infrared vari- 
ability of YSOs will provide key insights into the 
structure of YSOs and the physical mechanisms 
of the star formation process. 

In this paper, we present initial results and 
source lists of SAGE point sources for which we 
have detected thermal infrared variability between 
the epoch 1 and 2 photometry measurements. In 
§2, we outline the method of our approach to de- 
tect this variability. In §3, we discuss the results 
of the variable source selection and their identifi- 
cation of evolved stars and YSOs. In S4, we dis- 
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cuss implications of the results and summarize the 
main conclusions in §5. 

2. Method: Source Selection 

SAGE i s a Legacy project o n the Spitzer Space 
telescope ([Werner et al.l 120041 ) . which mapped a 
7° X 7° region of the Large Magellanic Cloud 



(LMC) using the IRAC camera in the [3.6 
[5.8], and [8.0] micron bands (jFazio et al 



[4.5' 



1200 



and the MIPS camera in the [24 ], [70], and [160] 



micron filters ( Rieke et al. I l2004l) . The scientific 



goals of the survey focus on the life cycle of bary- 
onic matter, as traced by dust emission, from its 
start in the interstellar medium (ISM), to the for- 
mation of new stars, to the death of these stars 
and the return of matter to the ISM. The sur- 
vey was done over two epochs with a total ob- 
serving time of 291 hrs with IRAC and 217 hrs 
wi th MIPS. The details o f the survey are described 
(|2006l ). 



Meixner et al 



The SAGE photomet- 
ric data have been separately extracted from the 
two epochs, that are separated by 3 months. The 
IRAC epoch 1 data were taken on July 15 - 26, 
2005, epoch 2 on Oct 26 - Nov 2, 2005 and the 
MIPS epoch 1 data were taken on July 27- Aug 3, 
2005 and epoch 2 on Nov 2 - 9, 2005. The SAGE 
Epoch 1 point source catalogfH containing over 
4 million sources of IRAC (SAGEcataloglRACe- 
pochl) and over 40,000 sources of MIPS 24 micron 
(SAGEcatalogMIPS24cpochl) have been merged 
with 2MASS JHK, and the Magellanic Clouds 
Photometric Survey (jZaritskv et al ]|2004|) . Analy- 
sis of the colors and magnitudes of this catalog has 
revealed three general types of LMC point sources: 
stars without dust, dusty evolved s tars and young 
stellar objects ( Meixner et al.|[2"006l ). Further clas- 
sificat ion of the SAGE sources by iBlum et"al 



( 20061 ) separated out the dusty evolved star classes 
such as supergiants, and asymptotic giant branch 
(AGB) stars from the red giants. The young stel- 
lar obje ct (YSO) population h as been identified in 
part bv lWhitnev et all (|2008[) . 

The SAGE 2-epoch point sources catalogs 
for IRAC (SAGEcataloglRACepochl and SAGE- 
catalogIRACepoch2) and MIPS (SAGEcatalog- 



^Catalogs that combine 2MASS JHK l lSkrutskie et al.l 
120061) . IRAC and MIPS 24 fim data are avail- 
able at the Spitzer Science Center website, 
http: / /ssc. spitzer. caltech.edu/legacy/all. html 



MIPS24epochl and SAGEcatalogMIPS24epoch2) 
were used for this study. Variable sources were 
identified using the matching and selection crite- 
ria described below and a variability index. We 
define the variability index V as the error-weighted 
flux difference in each SAGE band: 



V 



ifl - f2) 



where /i and /2 arc fluxes in the two epochs and 
cr/i and cr/2 are the associated errors. We also 
deflned a fractional flux (fF) for each band, the 
fractional change in flux 



1 1 - /2 
(/i + h)/2 



where /i and /2 are fluxes in the two epochs. The 
criteria for inclusion of sources into the catalog and 



the ar chive version are detailed in iMeixner et al 
(|2006l ). The source lists are stored in the SAGE 



database, a system implemented with Microsoft 
SQL server 2000. We developed scripts in Struc- 
tured Query Language (SQL) for creating the 
inter-epoch matching sources and for subsequent 
checks on, and quality control of, the variable 
source lists. 

2.1. Inter-epoch Matching of Sources 

AU four IRAC bands and the MIPS 24 micron, 
2-epoch data were used for this study. To find 
the variable IRAC sources in SAGE, the IRAC 
epoch 1 and epoch 2 catalogs were matched using 
a 2" search radius. Sources with multiple matches 
and sources with any neighboring sources within 
3" were excluded as careful consideration revealed 
that the PSF characteristics of the IRAC detectors 
make the measured fluxes of sources with neigh- 
bors closer than 2" less reliable. Furthermore, 
to avoid possibility of a mismatch between the 
epochs, only matches within 0.9" (the 3-(T value 
for IRAC photometry) were retained. We also 
used the following magnitude cuts to only include 
sources with highly reliable fluxes in our lists: 16, 
16, 14 and 13.5 mag at 3.6, 4.5, 5.8, and 8.0 ^m 
respectively. 

For MIPS, the 24/im PSCs SAGEcatalog- 
MIPS24epochl and SAGEcatalogMIPS24epoch2 
were used for this study. We used a 1" matching 
between the 2 epochs to eliminate wrong inter- 
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epoch matches. Both the individual epoch cata- 
logs have their astrometry matched to IRAC cat- 
alogs, which have been matched to the 0.3" accu- 
racy 2MASS catalog. We also note here that the 
MIPS24 catalogs are conservative catalogs that 
report only sources with reliable fluxes. With 
only two epochs of data we do not aim to pro- 
duce a complete list of infrared variable sources 
in the LMC; nonetheless we aim for a highly reli- 
able one, and require variability in more than one 
band. Figure [T] shows a histogram of the variabil- 
ity |V^| for the inter-epoch matched sources for all 
the IRAC and MIPS 24 jj,m bands. \V\ < 3 proba- 
bly indicates "random" errors in flux and |y| > 3 
a population being dominated by a systematic 
trend. Different bands are shown in different col- 
ors and the vertical lines at \V\ = 3 indicate the 
reliability cuts for the SAGE variables. 
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Fig. 1. — Histogram of all the SAGE sources 
matched between 2 epochs for the different bands 
of IRAC and MIPS 24 /Ltm band. Sources with 
\V\ > 3 in at least 2 consecutive bands are classi- 
fied as variable sources. 



2.2. Variability Criteria 

Our criterion for a variable source is > 3 in 
at least 2 consecutive bands in the same direction. 
The IRAC-selected sources were matched with 
the inter-epoch matched 24 fim sources. 24 fim 



sources with IV24I > 3 as well as IVs.ol > 3 in the 
same direction are deemed 24 fj,m variables. Ta- 
ble [1] explains the criteria used to identify the vari- 
able sources. The first part of the table character- 
izes the source catalogs. For each band the num- 
ber of sources with valid fluxes in the band at both 
epochs, the number and percentage of sources hav- 
ing IVbandl > 3, and the number and percentage of 
sources expected statistically to have iVbandl > 3 
are given. The second part describes the selec- 
tion criteria and the statistical significance of the 
variable sources. Here for each band, the sources 
not only have valid fluxes in both epochs for a 
given band but also in a neighboring band. Once 
again the numbers and percentages of sources both 
meeting the variability criterion and those that 
might do so based purely on statistical considera- 
tions (and therefore likely to be false variables) are 
given. Assuming that the measurement of flux for 
any source follows a gaussian distribution char- 
acterized by the error in the flux, the probabil- 
ity that two measurements of the flux of a non- 
variable source will not be within 3-sigma of each 
other is 0.27%. This applies for every band. We 
see that for all bands the number of sources with 
\V\ > 3 is more than that expected from statis- 
tical variations. Requiring that the sources have 
|y| > 3 in at least two bands in the same direc- 
tion considerably reduces the chances of random 
mis-classification as a variable source. Firstly a 
source has to have valid fluxes at both epochs 
in at least two neighboring bands, and secondly 
have \V\ > 3 in both of these bands in the same 
direction. Statistically this has a probability of 
0.0004% for bands with a single neighboring band 
and 0.0007% for bands with two neighbors. For ex- 
ample, at 24 fj,m less than 1 source would be clas- 
sified as a variable source due to random chance 
alone. The argument could be made that the flux 
measurement of sources does not follow a gaus- 
sian distribution, but this simple, idealized statis- 
tical approach gives us lower limits on the num- 
ber of false variables in the list. We perform other 
tests on the variable source list that do not rely on 
normal error distributions and uncorrelated errors, 
and those tests suggest that most of the sources 
in the list are really variables. A few sources have 
band-to-band inconsistencies in V, that indicates 
wrong photometry in one band. However, in these 
cases the the inconsistent band is not responsi- 
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ble for the variability criterion and these fluxes 
have been nulled in the tables. Of the ~ 4 million 
sources in the SAGE Catalog, Table [2] catalogs the 
properties (coordinates, fluxes and variability) of 
the 1,967 SAGE variables. 

3. Source Classification 

Using the variability criteria discussed in §2 we 
find 1,967 SAGE variables at IRAC and MIPS 
24 fim bands. 514 sources are variable {\V\ > 3) 
in all the five bands. Table [2] lists the properties 
of the 1,967 SAGE variables. 

Most of the SAGE variables can be classified as 
AGB stars, a small number as YSOs and a num- 
ber of sources that are presently unclassified but 
could be OB stars, RGB stars, post-AGB stars, 
PNe, background active galaxies or other classes 
of variables like cepheids, RR Lyrae stars or WR 
stars. We choose three color-magnitude diagrams 
(CMDs) to classify the SAGE variable sources. 
In Figure [2] all the sources in the SAGE epoch 
1 catalog are used to plot t he He ss diagram shown 
in grayscale. iBlum et al. (|2006l ) showed that the 
[J] - [3.6] vs. [3.6] CMD was most useful to see 
the separation of the 0-rich, C-rich and extreme- 
AGB population. Over-plotted are the epoch 1 
magnitudes of the SAGE variable sources classi- 
fied into 0-rich, C-rich and ex treme AGBs based 

( 20061 ) 



on the classification scheme of Cioni et al 



andjBl um et al. (2006). The features in the under- 
lying C MD (graysca l e) ar e labeled 'A', 'B', and 



C (see Blum et al. 2006: Nikolaev fc Weinberg 



120001 ): 'A' being the OB star locus at the bluest, 
faint end of the diagram. The first prominent fin- 
ger (labeled 'B') corresponds to young A-G super- 
giants (SGs). The finger 'C consists mainly of 
foreground dwarfs and giants. The next finger to 
the right in the figure represents late- type (mostly 



M) S Gs and luminous, 0-rich M stars (IBlum et al 
20061 ). The rest of the region above the tip of 
the RGB are divided into 0-rich, C-rich and ex- 
treme AGB zones. Keeping these divisions in mind 
lets us classify most of the variable sources as 
evolved stars. To identify the YSO population, 
we cross-correlate the variable sourc es with the 
SAGE- YSO fist from IWhitnev et al 



( 2008h and 



find 29 YSO candidates. Sources that are not clas- 
sified as evolved stars or YSO candidates are plot- 
ted in yellow and will be discussed further in § 13.11 
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Fig. 2. — Color magnitude diagrams (CMDs) of 
the sources detected as variables in the SAGE sur- 
vey. The SAGE variables over-plotted on a Hess 
diagram showing the distribution of all the stars in 
the SAGE Epoch 1 catalog. AGB stars are clas- 
sified into C-rich, 0-rich and extreme AGB fol- 
lowing the scheme followed by Cioni et al. 2006 
and Blum et al. 2006. YSO candidates are from 
cross-correlation with the SAGE- YSOs classified 
by Whitney et al (2008). SEDs of sources plotted 
with solid diamonds are shown in Figures [7ll9l &fTT] 

Meixner et al. (2006) showed that the [8.0] -[24] 
vs [8.0] CMD was most appropriate for separating 
the dusty objects based on mass loss. Figure [3] 
shows the [8.0] - [24] vs [8.0] CMD for the SAGE 
variables, with the same color-scheme used for the 
classification as in Figure [2] Figure |4] shows the 
[4.5] - [24] vs [24] CMD for the SAGE variables. 
This CMD is the one of the most useful ones to de- 
lineate the LMC stellar population from the back- 
ground population. We see that most of the un- 
classified sources are among the redder and fainter 
population in this CMD and are likely to be exter- 
nal galaxies. These sources are further discussed 
in the next section 5 13.11 
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Table 1 

Criteria for variable sources 



Criteria 


3.6 fim 


4.5 fim 


5.8 fiva 


8.0 fitn 


24 fim 


# of sources with valid fluxes 


2,962,938 


1,939,173 


354,290 


133,655 


25,412 


# with \Vband\ > 3 


15341 


9857 


2395 


1647 


1625 


% with \Vband\ > 3 


0.52 


0.51 


0.68 


1.23 


6.39 


# expected statistically 


8000 


5236 


957 


361 


69 


% expected statistically 


0.27 


0.27 


0.27 


0.27 


0.27 


# of sources with valid fluxes 












in neighboring bands 


1,787,038'' 


2,121,027*^ 


458,607" 


124,618'* 


13901° 


# meet variability criterion 


1520 


1827 


1582 


1316 


623 


% meet variability criterion 


0.09 


0.09 


0.34 


1.06 


4.48 


# expected statistically 


8 


15 


3 


<1 


<1 


% expected statistically 


0.0004 


0.0007 


0.0007 


0.0004 


0.0004 



''Valid fluxes at 3.6 and 4.5 nm 

''Valid fluxes at 3.6 and 4.5 /im or at 4.5 and 5.8 fim 
"^Valid fluxes at 4.5 and 5.8 /im or at 5.8 and 8.0 fim 

Valid fluxes at 5.8 and 8.0 fim 
"Valid fluxes at 8.0 and 24 fim 



3.1. Unclassified Sources 

Using 2MASS and 3.6 fim color criteria most of 
the variable sources are classified as evolved AGB 
stars (with 0-rich, C-rich and ext reme subclasses). 
Using the 'Whitnev et all (|2008l ) YSO candidate 
list we also identify YSO candidates as infrared 
variables. The classification of the variable sources 
that do not fall under any of these subcategories 
is more challenging. Approximately 17% of the 
variables remain unclassified and ~ 70% of them 
remain so mainly because they are missing IRAC 
and/or 2MASS fluxes. This could be because a) 
they are undetected being dust enshrouded, e.g 
very dusty YSOs, b) extremely bright and there- 
fore saturated in IRAC, or c) spatially extended 
at IRAC wavelengths and thus missing from the 
point source catalog. In preliminary follow-up 
work we have identified 20 unclassified sources as 
cepheids detected in the MACHO survey. Some 
of the other 'unclassified' variables fall on distinct 
regions of the CMD and could be variable OB 
stars, RGB stars or 0-rich stars that were too faint 
to be detected in the 2MASS survey. The other 
sources could be YSOs, Post-AGB sources, PNe, 
background galaxies or other classes of variables 



like cepheids, RRLyrae stars, WR stars, QSOs or 
AGN. Further follow-up work is be required to 
identify the specific nature of these sources. 

4. Discussion 

In this first paper we identify infrared variable 
sources in the LMC and provide an approximate 
classification as to their nature as discussed in §3. 
In this section, we discuss the nature of these 
sources and some implications of their infrared 
variability for SAGE studies. In particular, we fo- 
cus on the evolved AGB star population and YSO 
candidates. The total number of variables identi- 
fied as SAGE variables is 1,967 or ~0.06% of the 
~3 million IRAC point sources that are common 
to both epoch 1 and 2 catalogs . These fractional 
numbers suggest that we preferentially detect as 
variable the redder or dustier sources which tend 
to be AGB stars or YSOs. In Table [3] we sum- 
marize the different categories of sources we have 
detected as variables and compare them to other 
work done on t he classification o f SAGE sources as 
evolve d stars ( Blum et al. [ | 2006l) and YSO candi- 



dates (jWhitnev et al.ll2008l) and discuss them be- 
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Table 2 

Properties of the 1967 SAGE Variables 



IRAC 
Designation 


MIPS 
Designation 


E 


R.A. 

(dcg) 


DECL. 

(dcg) 


Class 




3.6 fj.ni 






4.5 fim 5.8 /iin 8.0 /jm 24 /iin 


flux 
(mjy) 


err 
(mJy) 


V 


fF 


... same eolumiis ... 


SSTISAGEIC ... 


SSTMISAGEI ... 


1 


79.9337 


-69.9940 


X 


108.80 


5.08 


4.1 


0.25 




SSTISAGE2C ... 


SSTM1SAGE2 ... 


2 


79.9335 


-69.9941 




85.05 


2.71 








SSTISAGEIC ... 


SSTMISAGEI ... 


1 


79.9421 


-68.8979 


Y 


0.72 


0.04 


-4.4 


-0.29 




SSTISAGE2C ... 


SSTM1SAGE2 ... 


2 


79.9422 


-68.8979 




0.96 


0.04 








SSTISAGEIC ... 




1 


79.9643 


-72.6359 


U 


4.93 


0.18 


-5.4 


-0.29 




SSTISAGE2C ... 




2 


79.9646 


-72.6359 




3.67 


0.15 








SSTISAGEIC ... 




1 


79.9718 


-69.6882 


o 


5.23 


0.18 


5.0 


0.29 




SSTISAGE2C ... 




2 


79.9719 


-69.6883 




6.97 


0.30 








SSTISAGEIC ... 




1 


79.9811 


-69.3925 


C 


25.27 


0.63 


3.3 


0.10 




SSTISAGE2C ... 




2 


79.9813 


-69.3925 




22.83 


0.39 









Note. — Table [2] is published in its entirety in the electronic edition of the Astronomical Journal. A portion is shown here for guidance regarding its form and content. 
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Fig. 3.— CMD of the SAGE variable sources. 
Classification of the AGB stars and underlying 
Hess diagram are same as that in the previous fig- 
ure. SEDs of sources plotted with solid diamonds 
are shown in Figures [7][n] & fTTl 



Unclassified So"urce's 
C-rich AG'B stars 
■ &-rich AGB stars t 
Extreme AGB stars c 
YSO condi.dates 




[4.5] - [24] 



Fig. 4. — CMD of the SAGE variable sources. 
Classification of the AGB stars and underlying 
Hess diagram are same as that in the previous fig- 
ure. SEDs of sources plotted with solid diamonds 
are shown in Figures [7][9] & fTTl 



low. Histograms (Figure [5]) of average variability 
over all bands show how the AGB stars, YSOs and 
unclassified sources differ in their variability. 



"T" 



Unclassified Sources 
C-rich AGB stars 
0-ricti AGB stors 
Extreme AGB stars 
YSO candidates 




20 



Fig. 5. — Histogram showing the variability distri- 
bution among the different classes of SAGE IRAC 
variables. Average variability in all 4 bands of 
IRAC and MIPS 24 ^m band has been plotted. 
C-rich stars in green, 0-rich stars in blue, extreme 
AGBs in red, YSO candidates in cyan and the un- 
classified sources in black. 

4.1. Evolved stars 

As expected, the variable source population is 
dominated by evolved stars, in particular the AGB 
stars. The SAGE variables have > 81% AGB 
stars, classified as either 0-rich, C-rich or extreme 
AGB stars. The spatial distribution of the SAGE 
variables, as shown in Figure [SI correlates approx- 
imately with the LMC 3.6 /im image which traces 
predo minantly the old stellar population of the 
LMC ( Blum et al. 2006f) and supports our iden- 
tification of the variables as evolved stars. In par- 
ticular, the highest density of sources is located in 
the main stellar bar of the LMC. The color cod- 
ing of the point sources on this image follows that 
of the CMDs. The extreme AGB stars, shown 
in red, dominate both the spatial distribution of 
variable sources and the CMDs (Figures [5] - H]) . 
Looking at the numbers of sources detected per 
category (Table [3]) we find that the extreme AGB 
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Table 3 
SAGE Variable Population 



Source Type 


SAGE Variables 


SAGE Epoch 1 Sources 


% detected as variable 


Total 


1967 


4,338,548 


0.05 


0-rich AGB 


353 


17,875^ 


1.975 


C-rich AGB 


426 


6,935^ 


6.14 


Extreme AGB 


820 


1,240'^ 


66.1 


YSO candidate 


29 


990*^ 


2.93 


Unclassified 


335 


4,311,508 


0.008 



^Classified by Blum ct al. (2006) 
b Classified by Whitney et al. (2008) 



stars outnumber the 0-rich or C-rich AGB stars 
by a factor of two. Moreover, the percentage of 
variables changes dramatically across the classified 
AGB sources. Only 2.0% of the 0-rich AGB stars 
have been detected as variable. This percentage 
increases to 6.1% for the C-rich AGB stars. For 
the extreme AGB stars, the percentage jumps to 
66%. This increase in percentage follows an evo- 
lutionary trend of the AGB. All stars on the early 
AGB start off as 0-rich. As they evolve to higher 
luminosities their pulsational periods increase, and 
many become C-rich AGB stars. At the highest 
luminosities, the periods are the longest and these 
extreme AGB stars are characterized by signifi- 
cant circumstellar dust emission. Thus, the frac- 
tional number of AGB stars identified as variable 
increases as the star becomes more evolved on the 
AGB. 

This preferential detection of the evolved AGB 
stars is similar to that in the IRAS variability 
studies in our Galaxy because our sampling period 
of 2 measurements separated by 3 months is simi- 
lar to the IRAS sampling. As iHarmon fc Gilmord 
(jl988l ) found in our Galaxy, we are preferentially 
finding the more evolved, dusty AGB stars which 
probably have longer periods. This preference is 
at least in part due to the sensitivities of detec- 
tion and of time sampling. Because they are more 
luminous and dust enshrouded, the extreme AGB 
stars are easier to detect and their variability is 
easier to measure reliably in the IRAC and MIPS 
bands. In addition, the period of the variability 
for the extreme AGB stars, ^400-500 days, is well 
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Fig. 6. — The spatial distribution of the variable 
sources in the LMC classified as AGB stars plotted 
on the 3.6 /xm map. Note the correlation of the 
spatial distribution of the variable population with 
the stellar density function (3.6 ^m). This con- 
firms that most of these sources are indeed stars. 
The colors represent the same classification of the 
sources as used in Figure [H 

matched to our two sample points in time sepa- 
rated by 3 months (^90 days) such that we could 
more likely than not measure a detectable fiux dif- 
ference in the two epochs. 
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Even though we do not detect all the AGB stars 
as infrared var iables, we do know that they are all 
variables (e.g. iFraser et al. 2005 ). The detected 
infrared variability reported here provides a con- 
straint on the possible range of the infrared flux 
in these types of AGB stars. The fractional differ- 
ence in fluxes we detect has a range of 0.1 to 1.5 
of the total flux. These changes in infrared flux 
will have significant implications for the analysis 
of SAGE data on the AGB stars. Estimates of the 
infrared excesses for this population will have some 
systematic errors because of the infrared variabil- 
ity. The mass-loss rates of the evolved stars will 
be derived from model fits to the spectral energy 
distributions of these sources, and the IRAC and 
MIPS 24 photometry provide constraints on 
the dusty shells. The infrared variability of the 
three different classes of AGB stars is revealed 
in the 2-epoch spectral energy distributions of a 
few 0-rich, C-rich and extreme AGB stars selected 
from our sample (Figures [7] - [H]) . The SEDs show 
the SAGE epochl data connected by dotted lines 
and epoch2 data by dashed lines. The 2MASS 
data are connected by solid lines and in many in- 
stances appears disjoint from either of the epochs, 
which is due to infrared variability as these data 
represent a different phase far removed from ei- 
ther of the SAGE epochs. The 0-rich and C-rich 
AGB stars show a similar amount of variation and 
the extreme AGB star shows a larger variation. 
Indeed the extreme AGB stars typically have a 
larger variability index than the 0-rich and C-rich 
AGB stars (Figure [5]). 

Some of the brightest variables are not classified 
because they are saturated in IRAC 3.6 /im band 
or 2MASS bands, which we use for classification 
purposes. The drop in fiux in the Raleigh- Jeans 
tails of the SED enables them to be detected at 
the longer IRAC bands (4.5, 5.8 and 8.0 /xm) and 
MIPS 24 fim without saturation. 

4.2. YSO candidates 

We have discovered a YSO candidate popu- 
lation with infrared variability. That this phe- 
nomenon was not noticed earlier is not surprising 
given that little work has been done on infrared 
variability of YSOs and most of it has been fol- 
lowup of known optically variable sources. These 
variable YSO candidates comprise at least 3% of 
all the YSO candidates. We say at least, be- 
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Fig. 7. — Spectral energy distributions showing 
the two epochs of a sample of variable 0-rich AGB 
stars. SAGE epoch 1 fluxes are connected by 
dotted lines and epoch 2 by dashed lines. The 
sources are marked with solid, blue diamonds on 
the CMDs. The 2MASS data are connected by 
solid lines and in many instances appears disjoint 
from either of the epochs, which is due to infrared 
variability as these data represent a different phase 
far removed from either of the SAGE epochs. 



cause there are a large number of sources in the 
same CMD space that are "unclassified" . The spa- 
tial distribution of YSO candidates correlates spa- 
tially with the LMC MIPS 24 micron emission, 
which traces massive star formation, supporting 
the identification of these variables as YSOs (Fig- 
ure [TO]). Also interesting is the spatial distribution 
of the "unclassified" sources (Figure [TOl) . These 
sources also closely follow the the 24 ^m emission. 
Further follow up work on extending the list of 
YSO variables that is in progress will help clarify 
if some of these unclassified sources are YSOs. 

In this study, the variability for this YSO popu- 
lation shows most prominently in the MIPS 24 /xm 
band (24 of 29 YSO candidates have 1^241 > 3). 
The reason for the predominance of MIPS 24 /im 
band variability becomes clear when we exam- 
ine the SEDs of a few variable YSO candidates 
(Figure ini). The SEDs show the SAGE epochl 
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Fig. 8. — Spectral energy distributions showing 
the two epochs of a sample of variable C-rich AGB 
stars. The sources are marked as solid, green di- 
amonds on the CMDs. SAGE epoch 1 fluxes are 
connected by dotted lines and epoch 2 by dashed 
lines. The sources are marked with solid, red dia- 
monds on the CMDs. The 2MASS data are con- 
nected by solid lines and in many instances ap- 
pears disjoint from either of the epochs, which is 
due to infrared variability as these data represent 
a different phase far removed from either of the 
SAGE epochs. 



data connected by dotted lines and epoch2 data 
by dashed lines. The 2MASS data are connected 
by solid lines and in many instances appears dis- 
joint from either of the epochs, which is due to 
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Fig. 9. — Spectral energy distributions showing 
the two epochs of a sample of variable extreme 
AGB stars. SAGE epoch 1 fluxes are connected 
by dotted lines and epoch 2 by dashed lines. The 
sources are marked with solid, red diamonds on 
the CMDs. The 2MASS data are connected by 
solid lines and in many instances appears disjoint 
from either of the epochs, which is due to infrared 
variability as these data represent a different phase 
far removed from either of the SAGE epochs. 



infrared variability as these data represent a dif- 
ferent phase far removed from either of the SAGE 
epochs. Most of these sources are very faint in the 
2MASS and IRAC bands, and increase rapidly in 
flux at the longer wavelengths. In fact, while our 
SEDs end at 24 /um the trend suggests that for 
some of these YSOs, the SEDs actually peak at 
even longer wavelengths. Thus, our sensitivity to 
these sources is greatest at 24 /im and the lower 
number of IRAC variable YSO candidates in the 
SAGE data may be due to the lower sensitivity in 
the IRAC bands to detect the variability. 

Whitnev et al. ( 2008h have completed a candi- 



date list of over 990 YSO candidates identified in 
the SAGE epoch 1 catalog. This list is a lower 
limit to the total number of YSOs in the LMC be- 
cause the selection was conservative and avoided 
regions in CMD space densely populated by other 
types of sources, e.g. AGB stars and background 
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Fig. 10. — The spatial distribution of variable 
massive YSOs and the unclassified sources plot- 
ted on SAGE 24 jim image. The YSO candidates 
are shown in red and the unclassified candidates 
in cyan. 

galaxies. Our variable YSO candidates comprise 
- 3% of the larger 990 YSO candidate Hst. It 
is difhcult to assess if this percentage is a lower 
or upper limit because the number of YSOs and 
variable YSOs is rather uncertain in different di- 
rections. Nevertheless, the percentage is signifi- 
cant enough to catch our attention as a new class 
of infrared variables. With only two epochs of ob- 
servation spread 3 months apart, we cannot dis- 
tinguish between a periodic or a more stochastic 
phenomenon for the YSOs. Clearly, further study 
of these sources is required to understand the na- 
ture of their variability. 

5. Conclusions 

We have detected infrared variable stars in the 
LMC using the SAGE survey. The variable source 
population is dominated by evolved stars, in par- 
ticular the AGB stars (81% of the SAGE vari- 
ables). We preferentially find the more evolved, 
dusty AGB stars which probably have longer pe- 
riods. We also present the discovery of a YSO can- 
didate population with infrared variability. These 
variable YSO candidates comprise about 3% of all 
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Fig. 11. — Spectral energy distributions (SEDs) 
of a sample of YSO variable sources. The sources 
are marked as solid, cyan diamonds on the CMDs. 
SAGE epoch 1 fluxes are connected by dotted lines 
and epoch 2 by dashed lines. The sources show 
SEDs consistent with embedded objects with dust 
emission. The sources are marked with solid, red 
diamonds on the CMDs. The 2MASS data are 
connected by solid lines and in many instances ap- 
pears disjoint from either of the epochs, which is 
due to infrared variability as these data represent 
a different phase far removed from either of the 
SAGE epochs. 

the SAGE variables. 
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